Cerenkov-like radiation in a binary Schrodinger flow past an obstacle 
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We consider the dynamics of two coupled miscible Bose-Einstein condensates, when an obstacle 
is dragged through them. The existence of two different speeds of sound provides the possibility for 
three dynamical regimes: when both components are subcritical, we do not observe nucleation of 
coherent structures; when both components are supercritical they both form dark solitons in one 
dimension (ID) and vortices or rotating vortex dipoles in two dimensions (2D); in the intermediate 
regime, we observe the nucleation of a structure in the form of a dark-antidark soliton in ID; the 
2D analog of such a structure, a vortex- lump, is also observed. 



Introduction. In the past few years, there has been an 
increasing number of studies of multi-component Hamil- 
tonian systems. This has been triggered primarily by 
the development of theoretical and experimental results 
in coupled atomic Bose-Einstein condensates (BECs) [T|, 
and of coupled nonlinear optical systems (where the 
coupling can be, e.g., between different polarizations 
of light or different frequencies) [2]. These have, in 
turn, motivated detailed mathematical investigations of 
such coupled systems, typically described by nonlinear 
Schrodinger (NLS) equations [3[. In the setting of BECs, 
that will be the primary focus of this study, mixtures 
of different spin states of 87 Rb [4] and 23 Na as well 
as two-component BECs with different atomic species, 
such as 41 K- 87 Rb [6] and 7 Li- 133 Cs @, have been cre- 
ated in experiments. In the same context, a wide variety 
of theoretical studies have examined ground-state solu- 
tions [8| and small-amplitude excitations [9], as well as 
the formation of other nonlinear structures such as do- 
main walls [TO], one-dimensional (ID) bound dark-dark 
and dark-bright soliton complexes , spatially periodic 
states 
waves 



12 1, vortex dipoles [13[, vortex rings and slaved 
14 1, coupled vortex lattices [10], an d so on. 



At the same time, many theoretical and experimental 
studies deal with the dragging of an "impurity" (e.g., a 
blue-detuned laser beam) through a one- component con- 
densate. This setting has been demonstrated to be pro- 
totypical for dark soliton formation in ID [TO, [T3l , and 
for vortex formation in 2D [TO]. These nonlinear waves 
can be thought of as a type of nonlinear Cerenkov radi- 
ation that is emitted, when the motion of the impurity 
is supercritical with respect to the local speed of sound 
of the BEC. Recently, a combined experimental and the- 
oretical study of the Cerenkov emission of phonons by 
a laser obstacle was reported [19] : in a different study 
[10] , it has been shown that in the case of large obstacles 
(and for a supersonic flow of the BEC), the Cerenkov 
cone transforms into a spatial shock wave consisting of a 
chain of dark solitons [20]. In fact, this setting has been 
particularly relevant for the study of the breakdown of 



superfluidity (and emergence of dissipation) and the asso- 
ciated Landau criterion J2l|. Indeed, earlier experiments 
[22[ have demonstrated the onset of dissipation, when a 
blue-detuned laser beam moves through the BEC with 
velocities above a threshold. We also note in passing 
that the appearance of similar effects (e.g., the backward- 
propagating Cerenkov radiation) in photonic crystals [23[ 
is yet another illustration of the interest in this research 
direction. 

In the present work we study the dragging of a S-\ike 
obstacle in a two- component superfluid flow. If the com- 
ponents are assumed to be immiscible, then they will 
tend to phase-separate and the problem reverts to its 
single-component version. For this reason, we consider 
the case of two miscible components, which is particularly 
interesting due to the existence of two distinct "speeds 
of sound". In this setting, we find two critical speeds 
< < . For v < , we show that the impurity 
propagates without emitting Cerenkov radiation in the 
form of nonlinear waves. For < < v, both com- 
ponents are supercritical and the impurity emits gray 
solitons (in ID) propagating downstream in both com- 
ponents. However, the most interesting regime is the 
intermediate one, where one of the components is super- 
critical, yet the other is subcritical, leading to the spon- 
taneous formation of dark-antidark solitary waves previ- 
ously predicted (in a stationary form) in [24J . We demon- 
strate that when the strength of the impurity tends to 
zero, the critical speeds tend to the corresponding speeds 
of sound, yet we show how they deviate from these val- 
ues for finite impurity strengths. We also consider the 
2D case, where we also obtain the analog of the dark- 
antidark state in the form of a vortex-lump wave. 

The paper is structured as follows. We first present the 
theoretical framework, and calculate the critical veloci- 
ties. We then numerically investigate the ID (both for 
untrapped and trapped BECs) and the 2D case. Finally, 
we summarize our findings and present our conclusions. 

Theoretical Setup. We consider the following coupled 
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NLS equations, describing a quasi- ID binary BEC Q: 
idt^j = i^-\d 2 x + X> fe W 2 + t/fj, (1) 

where i/jj (j = 1,2) are the mean- field wavefunctions, and 
V GX t = Vi + V2 is the external potential, assumed to be 
composed by a repulsive potential of a blue-detuned laser 
beam, Vi, and a trapping harmonic potential, V2, i.e., 

Vi = Ae-^- vt ^^\ V 2 = (2) 

where A, e and v are, respectively, the strength, width 
and velocity of the laser obstacle, and Q the harmonic 
trap strength. The nonlinearity coefficients are chosen to 
be #ii : #12 : #22 = 1.5 : 1 : 1.03. Notice that two of them 
have the ratios that are typical for 87 Rb [4|], while the 
third is tuned to a different value (so as to ensure misci- 
bility since the standard value of g\\ = 0.97 would lead to 
immiscibility ) . The tuning can be achieved by means of 
a Feshbach resonance [25]. Moreover, throughout the pa- 
per we use the following parameter values for our numer- 
ical computations: chemical potentials \i\ = 1.2, H2 = 1, 
obstacle width e = 0.5, and harmonic trap strength 
Q = 0.02. The results do not change qualitatively for 
other parameter values. 

The uniform solutions of Eqs. (pQ) satisfy 

|^(0) |2 = Mlff22 ~ M2ffl2 |^(Q) | 2 ^2^11 ~ Mlffl2 ^ 

where A = gii#22 — #12- Expressing Eqs. (PQ) in the 
traveling- wave frame (i.e., x — > x — vt) and linearizing 

around these uniform states according to ipj = ip^ 
we obtain the equations for the small amplitudes ijj^p , 

with j, k e {1, 2}, k ^ j and c 2 = gjj(ipj°^) 2 > Writing this 
system as one of four first-order equations (for = 
ipj 1 ^ = (d/dx)ip^ and V^)? namely, 
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where 6 = g\2W{ W2 • Then, the 4x4 matrix has eigen- 
values A 2 = c 2 -v 2 ±^c 4 - A(4 0) 4 0) ) 2 , with c 2 = (c 2 + 
c|)/2. For stability, we need the eigenvalues to be real, 
hence A 2 > implies that v < < where the crit- 
ical velocities, corresponding to the two distinct speeds 

of sound, are given by vi 1 ' 2 "' = c 2 ± \J c A — A^ip^ip^) 2 - 
For the parameters mentioned above, = 0.34393 and 



(2) 

Vc = 1.04796. We expect that superfluidity will break 
down when the speed v of the defect overcomes these 
speeds; in fact, as argued in [l6|, [2l|, the actual criti- 
cal point should be expected to be lower than the above 
Landau prediction. 
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FIG. 1: Top panel: "Stability" boundaries for static dark soli- 
ton pairs: the critical velocities for the first and the second 
component are shown as a function of the impurity strength 
A. Middle panels: Space-time evolution of the density con- 
tours for the two components in moving coordinate frame with 
velocity v = 0.2 (the speed of the impurity); clearly the im- 
purity induces the radiation of dark-antidark pairs. Bottom 
panels: similar to the middle panels but with v — 0.3, where 
both dark-antidark dipoles and dark-dark pairs are emitted. 
Impurity parameters: A = 0.9 and e = 0.5. 

Numerical Results. We now turn to the numerical in- 
vestigation of the above setting. In Fig. [TJ we test the 
theoretical prediction for the existence of two critical ve- 
locities for the dynamical evolution in the two compo- 
nents. The top panel of the figure shows a relevant "bi- 
furcation diagram" , where the dependence of the critical 
velocities on the "strength" A of the impurity is numer- 
ically evaluated. Note that as the strength of the impu- 
rity tends to zero (A — > 0) one recovers the numerical 
values for Vc ' stated above. The critical velocities are 
computed by finding (i) the speed above which, ap- 
parently, one component is supercritical emitting dark 
solitons, while the other is subcritical emitting antidark 
solitons (i.e., bright solitons on a finite background) that 

(2) 

"accompany" the dark ones; and (ii) the speed Vc above 
which both components nucleate dark solitons (see also 
bottom panel of the figure for a space-time evolution of 
the density contour plot for the two components for such 
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FIG. 2: Space-time contour plots of the components' density, 
when the obstacle's velocity is v = 0.2 (top panels) and v = 
0.3 (bottom panels) in the trapped case (Q = 0.02). Impurity 
parameters: A = 0.9 and e = 0.5. 



a supercritical — in both components — velocity v). It is 
noteworthy that such structures had been predicted in a 
steady form in [2J], but, to our knowledge, this is the first 
demonstration of their dynamical nucleation. A further 

(2) 

observation is worth making about the case of v > Vc . 
Note that, especially for early times, the impurity (sta- 
tioned at x = in the computations of Fig. [1] performed 
in the co moving frame) initially emits structures that ap- 
pear more like dark-antidark dipoles, i.e., dark-antidark 
pairs in one-component coupled with antidark-dark pairs 
in the other component. This is again the first manifes- 
tation of such structures (to the best of our knowledge); 
however, we will make a connection below to their 2D 
analog that has been previously proposed [261 ]. 

The above phenomenology also persists in the presence 
of an harmonic trap, which is a more realistic setting for 
magnetically confined BECs. This is clearly shown in 
Fig. [2j where all parameters are the same as in the corre- 
sponding plots of Fig. m but incorporating an harmonic 
trap of frequency Q = 0.02. 

We now turn to the 2D case where the second spatial 
derivatives in Eq. (pQ) are substituted by the Laplacian 
and the impurity potential is replaced by its 2D coun- 
terpart Vi(x,y) = (A/4)exp (-(x - vt) 2 /2e 2 ) (tanh(i/ + 
w/2) + l)(tanh(— y + w/2) + 1), modeling a light sheet 
of strength A, width e, and size w (see elongated vertical 
bar in panels a and c of Fig. [3j). In our 2D simulations 
we took A = 0.9, e = 0.5 (i.e. same parameters as for 
the ID case), and w = 10,15,20. Given the similarities 
of the trapped and untrapped case in the relevant phe- 
nomenology, we only show the latter here. In Fig. [3] we 
illustrate the two regimes leading to vortex nucleation 
(the trivial regime for subcritical velocity in both com- 
ponents is not shown here): (A) v = 0.235 is subcritical 
in the first component, but supercritical in the second, 
resulting for the latter in a vortex state which is coupled 
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FIG. 3: (Color online). Final snapshots after vortex nu- 
cleation for different velocities of a running impurity of size 
w = 10 (elongated vertical bar in panels a and c). The two de- 
picted cases correspond to: (A) v = 0.235 is subsonic for com- 
ponent ipi and supersonic for component ip2- (B) v = 0.350 
is supersonic for both components. The top (bottom) panel 
corresponds to component ipi (^2). Left panels (a,c) show 
the square modulus of the solution together with the mov- 
ing impurity. The right panels (b,d) show the vorticity cj 
in the rectangular area depicted in the left panel counter- 
parts. Red/blue (top/bottom spots) corresponds to regions 
with positive/ negative vorticity. 



to a lump (a 2D structure on a finite background) in the 
first component. Notice that the presence of vortex states 
is clearly illustrated in all the figures contained herein, by 
means of the contours of the vorticity uo = V x v s , where 
v s = (^*V^ - ^V^)/i\^\ 2 is the velocity field. Such 
structures have been reported previously for gn = #22 in 
[26( | . (B) v = 0.235 is supercritical for both components. 
This results in the formation of a dipole state which con- 
tains a vortex-lump pair, coupled to a lump- vortex pair, 
in a form similar to the stationary states reported in [l3[ . 

The last case (supersonic in both components) is exam- 
ined in further detail in Fig. 0] for different velocities and 
widths of the quasi- ID obstacle. In the figure, the actual 
spatio-temporal evolution of the vorticity is shown. This 
clearly reveals the presence of a vortex dipole between 
the two components; moreover, this robust type of state 
appears to be clearly rotating, as time evolves. Further- 
more, it can be noted that the wider the obstacle, the 
more complex the ensuing vortex patterns will be, with 
multiple vortex pairs being emitted. 

Conclusions. We have considered the nucleation of co- 
herent structures by a moving obstacle in two miscible 
BEC components. In one spatial dimension, we identi- 
fied three different regimes: one without nucleation; one 
involving the nucleation of dark-antidark solitons previ- 
ously predicted in stationary form in [24]; and one pro- 
ducing dark-dark soliton pairs, as well as dark-antidark 
dipoles. The critical points between these regimes were 
numerically obtained and, consistently with the corre- 
sponding single-component theory, were shown to ap- 
proach the Landau criterion for impurity strength tend- 
ing to zero; they were systematically lower than that as 
this strength increased. It was shown that similar behav- 
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FIG. 4: (Color online). Vortex nucleation by a running im- 
purity. The panels depict 3D contour plots of the vorticity 
cj(x,y,t) for the different velocities and impurity size combi- 
nations: left panel: w = 15 and v = 0.35 and right panel: 
w = 20 and v = 0.325. Both cases correspond to impurity ve- 
locities that are supersonic for both components. The red/blue 
isosurfaces correspond to positively /negatively charged vor- 
tices in the first component. We also superimpose the vortic- 
ity isocontours of the second component where green/magenta 
correspond to positive/ negative vorticity. Note how rotating 
vortex dipole pairs between the two components are formed 
(left-most red-green and blue-magenta intertwined vorticity 
lines) . 



ior also occurs in the case of the harmonically trapped 
coupled BECs. We also examined the same type of be- 
havior in 2D systems. We observed the existence of sim- 
ilar types of regimes, as in the ID case (subcritical in 
both, supercritical in one, and supercritical in both). 
The intermediate regime gave rise to vortex-lump type 
structures (also discussed in [261]), while the supercritical 
regime gave the first example of nucleation of vortex- 
lump dipoles (obtained in stationary form in [l3[), which 
were actually observed to be rotating as time evolved. 
This investigation indicates that there is an interesting 
spectrum of dynamical possibilities available in multi- 
component condensates, which it would certainly be rel- 
evant to explore experimentally. The recent realization 
of spinor condensates with more than two components 
may provide a fertile ground for further theoretical in- 
vestigations in such higher-component settings. 
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